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Abstract. We study the compositions of an integer n whose part sizes do not exceed a fixed 
integer k. We use the methods of analytic combinatorics to obtain precise asymptotic formulas 
for the number of such compositions, the total number of parts among all such compositions, 
the expected number of parts in such a composition, the total number of times a particular 
part size appears among all such compositions, and the expected multiplicity of a given part 
size in such a composition. Along the way we also obtain recurrences and generating functions 
for calculating several of these quantities. Our results also apply to questions about certain 
kinds of tilings and rhythm patterns. 



1. Introduction 

Let k be a fixed positive integer. We study numerical and combinatorial aspects of the 
sequence {F n }% ) = _ 00 , defined by 

{0, ifra<0; 
1, ifn = 0; (1.1) 

ELi F n-a, ifn>0, 

and related sequences. When k = 2, the numbers F n form the (shifted) Fibonacci sequence. 
The numbers F n are known to count many things, including the tilings of the 1 X n rectangle 
using tiles of length no more than [2j [3] , the rhythm patterns of length n with note lengths 
not exceeding k (Section [2j [7]), the new cells after the nth step in a budding process that 
starts with one cell and where each cell has one child per step, up to a total of k children per 
cell [15j, and the compositions of an integer n with part sizes no larger than k (Section[2j [8]). 

A composition of a positive integer n is a sequence of positive integers (Ai, A2, • • • , A m ) such 
that YZjLi^j = n - The f° ur compositions of n = 3, for instance, are (1,1,1), (1,2), (2,1), 
and (3). It is well known and easy to prove that a positive integer n has 2 n_1 compositions 
|10j . The integers Xj of a composition are its parts. If L is a set of positive integers, we can 
consider the set of compositions of n such that every part of each composition lies in L. Such 
compositions are called restricted compositions or L- compositions of n. 

The literature on unrestricted integer compositions is vast, and restricted compositions have 
also received a fair amount of attention. For example, {1, 2}-compositions were studied in [T], 
{1, /^-compositions were studied in [3], and L-compositions in general were studied in [8]. See 
also [U QUI HH E21 H3J HB] . In addition, both restricted and unrestricted integer compositions 
are examples of sequence constructions in symbolic combinatorics, so several very general 
theorems from this field apply to them [5]. 
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In this paper we fix k, let L = {1, 2, . . . , k}, and study the set of L-compositions of n, i.e., 
the set of compositions of n whose part sizes do not exceed k. We use the following notation. 

• F n denotes the number of {1, 2, . . . , /c}-compositions of n. (We show that this definition 
of F n agrees with (jl.ip in Section [2j) 

• T n denotes the total number of parts among all {1, 2, . . . , /c}-compositions of n. 

• A n = T n /F n denotes the average number of parts of a {1, 2, . . . , /c}-composition of n. 

• For j G {1,2,..., k}, C n j denotes the total number of times the part size j appears 
among all {1, 2, ... , fc}-compositions of n. 

• For j £ {1, 2, . . . , k}, A n j = C n j/F n denotes the average number of times the part 
size j appears in a {1, 2, . . . , /cj-composition of n. 

In this paper we study these five quantities in a unified way, with the ultimate goal of giving 
precise asymptotic formulas for all of them. When we say that our asymptotic formulas are 
precise, we mean that both their percentage errors and their absolute errors go to as n — > oo. 
In fact, the absolute errors for our asymptotic formulas decay at an exponential rate. 

Previously, Flores [6] gave a precise asymptotic for F n , which he derived using the theory 
of recurrence relations whose characteristic polynomials have no repeated roots. Then, a 
precise asymptotic for A n was obtained in [5, Theorem V.l] using the methods of analytic 
combinatorics. Combining these immediately results in a precise asymptotic for T n . Our 
asymptotic for C n j is new, and our asymptotic for A n j improves on the one in [5j Theorem 
V.2]. 

Along the way to developing our asymptotic formulas, we also give recursive methods for 
computing all five of these quantities and generating functions for F n , T n , and C n j, some of 
which have appeared previously. Since one of our goals is to provide a unified treatment of 
these five quantities, we provide proofs of all of our results, including the ones that are not 
new. In addition, our approach provides an alternate route to Flores's asymptotic for F n . 

We use the following notation. Let G(x) = Yla=i %a an d let <P denote the unique positive 
real solution to 

^V + --- + ^ = L M 
(That there is a unique positive real solution is guaranteed by Descartes' rule of signs.) Let G 1 
and G" denote the ordinary first and second derivatives of G, and let a = l/4>- When k>2 
we have 4> > 1, and when k = 2, 4> is the golden ratio. Here are the asymptotic formulas we 
derive in this paper. (We review big-O notation in Section |U) 

Theorem 1.1. For some constant A with < A < 1 we have 

a) p n = |^ + o{a») n> 

(3) An = G^) {n + 1} " 1 + + ° {An) The ° rem V - 1] ' 

(4) C n>j = ^(n + * G , {a)3 ( } + 0{A"), and 

(5) ^ = ^(„ + l- i) + ^) +0 (A-). 

Actually, the formula for A n in Theorem 11.11 follows from an asymptotic formula for su- 
percritical sequences, of which {1,2,... , fc}-compositions are one example. See [5] for more 
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on supercritical sequences. Previously, the best asymptotic result we could find for A n j is 
an asymptotic for supercritical sequences (5[ Theorem V.2] which, when translated into the 
setting of {1, 2, ... , fe}-compositions, becomes 

^nj = ^yn + 0(l). 

Our asymptotic formula for A n j refines this into a precise asymptotic. Also, we remark that 
the formula for F n in Theorem 11.11 was stated differently in [6] , but one may use equation (jl.2p 
to verify that the formula here is the same as the one in [6]. 

The rest of the paper is organized as follows. In Section [2] we review the connection between 
rhythm patterns and restricted compositions and we derive our recurrences. In Section [3] we 
use our recurrences to obtain our generating functions. We then analyze those generating 
functions in Section [5] to obtain our asymptotic results. Finally, Section [5] contains tables of 
values of F n , T n , A n , C n j, and A n j and our asymptotic approximations to them for small 
values of n and k. 

2. Rhythm patterns and recurrences 

In this section we give recursive methods for calculating the quantities F n , T n , A n , C n j, 
and A n j. Before doing so, however, we mention that our original motivation for studying 
restricted compositions was musical. The L-compositions of an integer n are in bijection with 
the rhythm patterns of length n with note lengths in L [7]. A rhythm pattern of length n is 
a sequence of notes (also called hits) and rests played over n evenly spaced pulses, where a 
note occurs on the first pulse, and on each further pulse either a note or a rest (but not both) 
occurs. The length of a note is 1 plus the number of rests until the next note or the end of 
the pattern. The bijection is given by mapping compositions to rhythm patterns part by part, 
where the part j corresponds to a note of length j (i.e., a hit followed by j — 1 rests). For 
example, the composition (2, 1, 1, 1, 2, 2, 1) of n = 10 corresponds to the rhythm pattern 

hit - rest - hit - hit - hit - hit - rest - hit - rest - hit. 

Thus, if we consider {1, 2, . . . , /c}-compositions, F n counts the number of rhythm patterns of 
length n with no note length exceeding k, T n counts how many notes one would play if one 
played all such patterns one after another, and A n says how many notes occur on average 
in such a pattern. C n j indicates how many times the note length j occurs among all such 
patterns, and A n j indicates how many times the note length j occurs on average in such a 
pattern. 

We now derive our recurrences. Our recurrence for F n is well known, and previously ap- 
peared in [51 [15], among others. 

Theorem 2.1. F n satisfies 

'0, ifn<0; 
1, ifn = 0; (2.1) 

Xtl^n-a, ifn>0. 

Proof. Certainly F n = if n is negative. Fq = 1 as there is one composition — the empty 
composition — of 0. If n > then we obtain the {1, 2, . . . , fc}-compositions of n by appending 
l's to the {1, 2, ... , fc}-compositions of n — 1, by appending 2's to the {1, 2, ... , /c}-compositions 
of n — 2, etc., all the way down to appending /c's to the {1,2,... , fc}-compositions of n — k. 
Thus we have F n = F n _i + F„_2 + • • • + F n _k for n > 0. □ 
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Next we obtain a recurrence for C n j (and hence a recursive method for computing A n j). 
Theorem 2.2. Let j € {1, 2, . . . , k}. Then 



Proof. If n < j we have C n j = 0, with which this recurrence agrees. Next, if n = j, since 
Fq = 1 and C a ,j = for a < j, this recurrence correctly yields C n j = 1. 

Finally, suppose n > j. By induction we may assume that this recurrence correctly computes 
C a j for all a < n. Split the set of {1,2,..., /c}-compositions of n into two classes: those that 
end with j and those that do not. Among those that end with j, there are (C n -jj + F n _j) j's. 
(Specifically, there are a total of F n _j j's at the ends and a total of C n -jj j's among all other 
places of these compositions.) Among the {1, 2, . . . , /c}-compositions of n that end with a € 
{1,2,..., k}, for a ^ j, there are (C n _ a j) j's. Summing these counts across a G {1, 2, . . . , k} 
yields C nj - = F n _j + J2a=i C n -a,j, as desired. □ 

Finally we obtain a recurrence for T n (and hence a recursive method for computing A n ^). 

Theorem 2.3. T n satisfies 



Proof. Tq = because there are no parts in the empty composition. If n > 0, then we have 




(2.2) 




(2.3) 



k 



T n — y c, 



0=1 




k 



k 




a=l i=l 



A- 




o=l 



using the fact that the recurrences C n j = F n _j + Ylt=i Cn-a,j and F n = Y2j=i ^ 
valid for n > 0. 



are both 
□ 



3. Generating functions 



In this section wg derive generating functions for F n , T n ^ and C n j. Recall = 
The generating function for F n is well known, and previously appeared in [5j[8]. 



o=l 



x 



,a 



Theorem 3.1. Let F(x) be the ordinary generating function for the number of {1, 2, . . . , A;}- 
compositions of n. That is, F{x) = Fq + F\x + F2X 2 + • • • . Then 

„, s 1 
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Proof. We prove this directly. Using ^0 = 1 and the recurrence in (I2.1D for n > 0, we have 

oo 

F{x) = l + Y J FnX n 



n=l 

oo k 



1 + 

n=l a=l 

k oo 

a=l n=l 
k oo 
1 + J> a X>n-aZ n - a 



a=l n=l 

fc 

= l + £VF(x). 

a=l 

Solving for F(x) yields F{x) = 1/(1 — X) a =l as desired. □ 

Next, for j € {1, 2, . . . , k}, we obtain the generating function for C n j. We remark that this 
generating function can also be obtained by differentiating formula (7) on p. 293 of [5] with 
respect to u and evaluating at u = 1. There the appropriate generating function was built 
directly from the combinatorics of sequences. Here we derive it from our recursive formula 
(1231) instead. 



Theorem 3.2. Let j £ {1, 2, . . . , k} and let Cj(x) be the ordinary generating function for the 
total number of occurrences of the part size j among all {1, 2, . . . , k} -compositions of n. That 
is, Cj(x) = Coj + C\jx + C2JX 2 + • • • . Then 

° 3{x) = (T^RT 

Proof. Using Coj = and equation ()2.2[) we have 

00 

Cj(x) = y c n jx n 

n=l 

00 / k \ 



^ ( Fn-j + ^2 C n -a,j 
71=1 \ a=l J 

00 00 k 

n=l n=l a=l 

00 fc 

X"' ^ F n —jX 3 + ^ ^ X ^ ^ C n — a jX 

n=l a=l n=l 
fe 



a=l 



Solving for Cj(x) we obtain Cj{x) = x 3 ' F(x)/{1 — Yla=i x<1 ) = x 3 /{\ — Y2a=l x )i as desired. 
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Theorem 13.21 also reveals the nice fact that for a fixed value of k, the sequences {C ni j}^L 
for j 6 {1, 2, . . . , k} are all just shifts of the same sequence. See Table [TJ for example. 

Finally we obtain the generating function for T n . This generating function previously ap- 
peared more generally in the context of sequence constructions on p. 178 of [5]. Here we derive 
it using our recursive formula (12.3|) instead. 



Theorem 3.3. Let T{x) be the ordinary generating function for the total number of parts 
among all {1, 2, . . . , k} -compositions of n. That is, T(x) = Tq + T\x + T^x 1 + • • • . Then 

{ 1 (l-G(x)) 2 l-G(x)' 
Proof. Using To = 0, the recurrence in (|2.3p . and Theorem 13 .!( we have 

T(x) -- 




1 - G(x) 

v ' a=l 

Solving for T(x) yields T(x) = 1/(1 - G{x)) 2 - 1/(1 - G(x)), as desired. □ 

4. ASYMPTOTICS 

In this section we use the analytic properties of our generating functions to obtain our 
asymptotic formulas for F n , T n , A n , C n j, and A n j. In this section we use n to stand for 
a nonnegative integer, x a real number, and z a complex number. In addition, for purposes 
of nondegeneracy in our arguments, in this section we assume k > 1. However, it is easily 
verified by inspection that all of our asymptotic results also hold for k = 1. In Section [4. II we 
review the ideas from analytic combinatorics and establish several facts we need to derive our 
asymptotic results. In Section I4T21 we prove our asymptotic results. 

4.1. Preliminaries. We are going to obtain our asymptotic formulas using the methods of 
analytic combinatorics [5j [17]. Namely, we will view the generating functions obtained in 
Section [3] as functions of a complex variable and we will analyze their dominant singularities 
to give us information about the growth rate of their coefficients. First, we quickly review 
big-O notation. 

Definition 4.1. Let Let N = {0, 1,2,.. .} and let f,g : N -> R. We write f(n) = 0(g(n)) if 
there are positive real numbers C and N for which 

|/(n)| < C\g(n)\ 

whenever n > N . 
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If f,g, and h are functions, whenever we write an equation like 

f(n) = h(n) + 0(g(n)), 

what we mean is 

f(n) = h(n) + E{n) 

for some function E, where E{n) = 0(g(n)). 

The basic observation from complex analysis that we will need is this |17t Theorem 2.4.3]. 

Theorem 4.2. Let f(z) = Y^ri=o a n zn ^ e analytic in a region containing the origin, and let 
zq be any singularity of f(z) of smallest modulus. Let R = \zq\ and fix a value e > 0. 
Then there exists N such that, for all n > N , we have 

k'J < I 4 + e 



Also, for infinitely many n we have 



R 



R 



The main idea is then the following (see also [17]). Let H(z) = ho + h\z + /i2^ 2 + • ■ ■ 
be a generating function (in particular, hi > for all i) with only isolated singularities and 
whose coefficient growth rate we'd like to understand. We use the standard notation [z n ]H(z) 
to stand for the coefficient h n of z n in H(z). Suppose that, as a function of a complex variable 
z, H(z) is analytic in a region containing the origin and has radius of convergence R. Then 
H(z) has a singularity on the circle {z S C : \z\ = R}. Suppose there is only one singularity zq 
on this circle (as there will be for our generating functions), and let S(z) be the principal part 
of H(z) at that singularity. We call zq the dominant singularity of H(z). Then the function 
H(z) — S(z) is analytic in a disk of radius R' > R centered at the origin. Thus, for any fixed 
e > 0, the coefficients of its expansion at the origin cannot grow faster than 



for sufficiently large n, and hence 



[ z ^H(z) = [z n }S(z) + 0^ + ey^ 



This is especially good when R' > 1: In this case we obtain that [z n ]S'(z) is a precise asymptotic 
for h n . (In particular, we obtain that [z n ]H(z) = [z n ]S(z) + 0(A n ) for some number A with 
< A < 1.) We remark that even when it is not the case that R' > 1, one can still use this 
formulation to obtain asymptotics — for many examples, see [5^ \17\. 

Recall that G(z) = Y^ a =i z<1 ' $ denotes the unique positive real solution to 

11 1 _ 

(ft 1 (ft 2 tfi k 

and a = l/4>. Thus G(a) = 1, so a is a singularity of F{z) and hence of all of our generating 
functions. In fact a is the dominant singularity of all of our generating functions, and we also 
have R' > 1 for all of our generating functions: 

Lemma 4.3. The only number z £ C for which G(z) = 1 and \z\ < 1 is z = a. 
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Proof. Since we have assumed k > 2, we have eft > 1, so < <r < 1. Miles [H] showed that 

M(z) = z k — — z — 1 has one root of modulus greater than one, and k — 1 (distinct) 

roots of modulus strictly smaller than one. It is straightforward to verify that r £ C is a root 
of M(z) if and only if 1/r is a root of Gr(z) — 1, which in turn yields the stated result. □ 

Thus, the expansion of our generating function F(z) at z = a will play a role in all of our 
asymptotic derivations. 

Lemma 4.4. The Laurent expansion of 

1 



1 - G{z) 

at z = a is 



F(z)- 
1 G"(a) 



G'(cr)(z — cr) 2G'(<7j^ 
/or some (unimportant) complex coefficients h,l2, ■ ■ ■■ 

Proof. Descartes' rule of signs implies that z = a has multiplicity one as a root of 1 — G{z) 
Thus the desired expansion is 

a . , , , s . , , , 2 



+ b + h(z - a) + l 2 (z - a) z + 



z — 0" 

for some constants a, b, li, We need to establish the values of a and 6. Write 1 — G(z) = 

(z — a)g(z) for some polynomial g(z) of which a is not a root. First, a is the residue of 
1/(1 — at z = a, so a = l/g(a). Note that, since a is a positive real number, the 

coefficients of g(z) are all real numbers. To compute a, then, we use the continuity of g along 
the real line and l'Hospital's rule to obtain 

, s ,. /x ,. g(x)(x - a) ,. 1 — G(x) -G'ix) . 
5(0-) = lim #(x) = hm - = lim — = hm — = -G'(a), (4.1) 

a;— ><7 x— ¥<j x — 0~ x— >-cr x — O x^ra \ 

and hence a = —1/G'(a). 
Next, we have that 

1 (z — a) (z — a) . , . .n , . .0 

- a + 6(z - o-) + /i(z - a) 2 + l 2 (z - a f + ■ ■ ■ 



g(z) (z-a)g(z) 1 - G(z) 
is analytic in a disk centered at z = a, so 

d ( 1 



dz \g(z) 



9>) 
5(0") 



2 ' 



We note that differentiating the equation g{x){x — a) = 1 — G{x) and performing a little 
algebra yields 

g'(x){x - a) = -G'{x) - g(x). 
To compute </(o~) we use the continuity of g' along the real line and l'Hospital's rule again: 

g'(x)(x-o-) - G '(s)-g(s) _ G "(x)-</(x) 
g (a) = hm g (x) = hm = hm = lim 



x— X — 0~ x—xt X — a x^ra 1 

= -G"(a)-g'(a). 

(The use of l'Hospital's rule here is justified by (|4.1|) . which says that lim^o- —G'{x) — g(x) 
—G'(a) — g(cr) = 0.) Finally, solving this equation for g'(cr) yields 
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and hence 

,_ G"{a) 



2G'{a) 2 ' 

as claimed. □ 

Finally, we will also need the straightforward facts that, expanding at the origin, we have 

1 (\ z z 2 \ , . 

Z + Z5 + ZS + -" > ( 4 - 2 ) 



i — (j \ r> a 2 a 3 
and 

1 1 2z 3z 2 , 

(z — a) z a z a d cr 4 

4.2. Asymptotic formulas. We are now ready to establish our asymptotic formulas. We 
begin with F n . The asymptotic for F n stated here is equivalent to the one in [6], where it 
was derived using the theory of recurrence relations whose characteristic polynomials have no 
repeated roots. Here we offer an alternate derivation. 

Theorem 4.5. We have 

<b n+l 

F ^hp) +0{An) 

for some number A with < A < 1. 

Proof. By Lemma 14.4} the expansion of the generating function F(z) = 1/(1 — G(z)) at its 
dominant singularity, z = cr, is 

~ nn J r + J(z-a) 

G'(a){z-a) 

where J(z — a) consists of (z — a) and higher-order terms. The expansion of the principal 
part of this function, 

-1 



G'(a)(z-aY 

at the origin is, by Equation (14.21) . 

1 f l z z2 

Thus 

[z n ]S(z) 



and so 

1 „ // 1 



for any e > 0, where R' is the distance from the origin to the singularity of F(z), other than 
a, closest to the origin. Lemma 14.31 implies that R' > 1, so e > can be chosen for which 
< (1/R' + e) < 1. Since <p = V* 7 ' this establishes our result with A = 1/R' + e. □ 

Next we turn to T n . 
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Theorem 4.6. We have 

& n+2 , , <p n+l G"(a) (j) n+1 

/or some number A with < j4 < 1 . 

Since A n = T n /F n , if one accepts as given the asymptotic for A n from Theorem 11.11 one 
can immediately deduce this asymptotic for T n by multiplying the asymptotic formula for F n 
in Theorem 14.51 by the one for A n in Theorem 11.11 We will also give our own proof of the 
asymptotic for A n below which is based on our asymptotic for T n , so we also offer the following 
direct derivation of Theorem 14.61 



Proof of Theorem \4-6\ By Lemma 14.41 we have that the expansion of the generating function 
T(z) = 1/(1 — G(z)) 2 — 1/(1 — G(z)) at its dominant singularity, z = a, is 

1 G'{af-G"{a) , 



G'(a) 2 (z-a) 2 G'{af(z- a) 

where J(z — a) consists of (z — a) and higher-order terms. The z n term of the expansion of 
the principal part of this function at the origin is, by Equations (|4.2p and (|4.3f) . 

(n + 1) + G "{a)-G'(a) 2 



G'(a) 2 a n + 2 G'iafa^ 1 ' 

so 

T - + | g »- g » 2 off 1 u 
n G'(a) 2 a n + 2 G'(a) 3 a n + l \\R' 

for any e > 0, where R' is the distance from the origin to the singularity of 1/(1 — G(z)), other 
than a, closest to the origin. Again, Lemma 14.31 implies that R' > 1, so e > can be chosen 
for which < (1/R' + e) < 1, establishing our result with A = 1/R' + e. □ 

Next we turn to A n , whose asymptotic here was proved in a more general setting in [5j 
Theorem V.l]. Our proof is similar to theirs and foreshadows our proof of Theorem 14.91 

Theorem 4.7. We have 

G'(o-y ' * ' G'(a) 
for some number A with < A < 1. 

Proof. From Theorems 14.51 and 14.61 we have 

Fn= G^) +E{n) = G^) { 1 + 
and 

r " = GK ( " + 1)+ G'(af -G^) +E2(n) 
for some functions E and with 2£(n) = 0(A n ) and E^n) = 0(A n ) for some constant ^4 
with < A < 1. Let I? be any constant with A < i? < 1. Then 

2; / G» _ / 1 \ 

n F n VG'(a) 1 ^ G'(a) 2 <^»+ 1 / V 1 + £(n)G'((7)/f + 1 J 

^<» +1 > + f$- 1+ H(i+iiw) 
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A n = ( 7^(" + !) + Tv7~k ~ 1 + °(^) ) (1 + °(^)) 



for 

_ £ 2 (n)GV) 


and 

Ea(n) _ W), 

since o~ is constant and > 1, it is immediate that E^{n) = 0{A n ) and E± = 0(A n ). It is 
therefore also straightforward to obtain 

—L— = l + 0(A n ). 
1 + E±(n) 

Thus, we have 

G"(a) 

G'(a) v " ' ^ ' G"(o-) s 

Multiplying this out and using nO(A n ) = 0(B n ) along with the fact that G'(a) > 0, G"(a) > 
0, and 4> > are constants, we obtain 

• 4 « = G^) (n+1) - 1+ i^ +0(B " ) 

Since < B < 1, we are done. □ 

Next we turn to C n j. 
Theorem 4.8. Let j e {1,2,... ,&}. T/ien we /jaue 

c -W , " +H1+1 ^r +0(i * ) 

/or some number A with < vl < 1. 
Proof. We have 

T i . I 

C n , = r*"l = f^-Jl 



lJ (1 — G(x)) 2 J (l-G(x)) 2 ' 

Lemma 14.41 implies that the expansion of 1/(1 — G{z)) 2 at its dominant singularity, z = a, is 

1 G"(a) 

G'(a) 2 (z-a) 2 G'(a) 3 (z-a) + {Z ah 

where J(z — a) consists of (z — o~)° and higher-order terms. The z n term of the expansion of 
the principal part of this function at the origin is, by Equations (14.21) and (|4.3j) . 

(n + 1) + G"(a) 



G'(a) 2 a n + 2 G'{afa n + 1} 

and thus 

C {n + 1 ~ j) I G " {a) I C(( l I e 

n ' 3 ~ G'(a) 2 o- n + 2 -J + G>(aYa^-i + \\R> + 

for any e > 0, where R' is the distance from the origin to the singularity of 1/(1 — G(z)), other 
than a, closest to the origin. Again, Lemma 14.31 implies that R' > 1, so e > can be chosen 
for which < (1/R' + e) < 1, establishing our result with A = 1/R 1 + e. □ 



Finally, W6 handle A n 



3- 
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Theorem 4.9. Let j e {1,2,... ,k}. Then 

A - 4-j (n I 1 i) I I Q(A») 

for some number A with < A < 1 . 

Proof. Our proof is similar to the proof of Theorem 14.71 By Theorems 14.51 and 14.81 we have 

F - <t> n+l . F(n) - ^ n+1 (, . 



and 

6 n+2 -i 6 n+1 -iG"(a, 



for some functions i£ and E<i with = 0(^4") and Ei(u) = 0(A n ) for some constant A 

with < A < 1. Let B be any constant with A < B < 1. Then 



F n VGV) G"(cr) 2 J \l + E(n)G' (a) /(/>"■ 



PS 



.(n+l-j) + Z G > } +g 3 (n 



for 



and 



£? 3 (n) 



An+1 



£(n)G'( ff ) 



0n+l ' 

since cr is constant and > 1, it is immediate that Ez(n) = 0{A n ) and E± = 0(A n ). It is 
therefore also straightforward to obtain 

' \ + 0{A n ). 



1 + E 4 (n) 
Thus, we have 

^ = (|^(n + 1 - j) + + 0(A»)J (1 + 0(A»)) . 

Multiplying this out and using nO(A n ) = 0(5") along with the fact that G'(a) > 0, G"(cr) > 
0, and 4> > are constants, we obtain 



A,,^(n + 1 -,) + ^ + (B.). 
Since < B < 1, we are done. 



□ 
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Table 1. Values of C n j for k = 2,3,4 
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k = 2 


n 


Cn,l 


Cn,2 











1 


1 





2 


2 


1 


3 


5 


2 


4 


10 


5 


5 


20 


10 


6 


38 


20 


7 


71 


38 


8 


130 


71 


9 


235 


130 


10 


420 


235 


11 


744 


420 


12 


1308 


744 


13 


2285 


1308 


14 


3970 


2285 


15 


6865 


3970 


16 


11822 


6865 


17 


20284 


11822 


18 


34690 


20284 


19 


59155 


34690 


20 


100610 


59155 



k = 3 


n 


C n i 


Cn.2 


C n 3 














1 


1 








2 


2 


1 





3 


5 


2 


1 


4 


12 


5 


2 


5 


26 


12 


5 


6 


56 


26 


12 


7 


118 


56 


26 


8 


244 


118 


56 


9 


499 


244 


118 


10 


1010 


499 


244 


11 


2027 


1010 


499 


12 


4040 


2027 


1010 


13 


8004 


4040 


2027 


14 


15776 


8004 


4040 


15 


30956 


15776 


8004 


16 


60504 


30956 


15776 


17 


117845 


60504 


30956 


18 


228818 


117845 


60504 


19 


443057 


228818 


117845 


20 


855732 


443057 


228818 



k = A 


n 


Cn,l 


Cn,2 




Cn 4 

















1 


1 











2 


2 


1 








3 


5 


2 


1 





4 


12 


5 


2 


1 


5 


28 


12 


5 


2 


6 


62 


28 


12 


5 


7 


136 


62 


28 


12 


8 


294 


136 


62 


28 


9 


628 


294 


136 


62 


10 


1328 


628 


294 


136 


11 


2787 


1328 


628 


294 


12 


5810 


2787 


1328 


628 


13 


12043 


5810 


2787 


1328 


14 


24840 


12043 


5810 


2787 


15 


51016 


24840 


12043 


5810 


16 


104380 


51016 


24840 


12043 


17 


212848 


104380 


51016 


24840 


18 


432732 


212848 


104380 


51016 


19 


877400 


432732 


212848 


104380 


20 


1774672 


877400 


432732 


212848 



5. Tables for small values of n and k 

In this section we give tables of values for F n , T n , A n , C n j, and A n j and our asymptotic 
approximations to them for small values of n and k. 

First, Table [TJ gives the values of C n j for j G {1,2,..., k}, for k = 2, 3, 4 and n = to 
n = 20. As we mentioned after the proof of Theorem 13. 2\ for a fixed value of k, the sequences 
{C n j}^L for j £ {1,2, ... , k} are all just shifts of the same sequence. 

Next, Tables[2l[3l andHJcontain the values of F n , T n , A n , C n> i, and A n ^\ and our asymptotic 
formula approximations to them for k = 2,3,4 and n = 0tora=15. Table [2] corresponds to 
k = 2, Table [3] corresponds to k = 3, and Table H] corresponds to k = 4. Entries in these tables 
are rounded to 3 decimal places. In these tables, the approximating formulas we are using are 
those given in Theorem 11.11 without the 0{A n ) terms. So, for instance, our approximation to 
F n is given by 

p ~ z. 

n G'(aY 

We also remark that for k = 2, Binet's formula leads to the same approximating asymptotic 
for F n — for n > 0, Binet's formula for the nth Fibonacci number yields 
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Table 2. Values of F n , T n , A n , C njl , and A nA for k = 2; 4> = (1 + \/5)/2 



71 


Tp 

r n 


Appr. 


± n 


Appr. 


A 

sin 


Appr. 


ri 

W,l 


Appr. 


A 

^in,l 


Appr. 


U 


1 


n 70 A 

u. r z4 





(J. (Joy 


n n 
0.0 


n i o a 
U.lz4 


n 
U 


U.l f\) 


n n 
0.0 


n o a v 


1 


1 


1 1 71 

1.171 


1 


n nno 


1 n 
1.0 


r\ o at 

(J. 847 


1 


n o i o 


1.0 


n en/1 


o 
Z 


o 

Z 


1 QCi A 


3 


z.y /b 


1.0 


1.0 ( 1 


o 
z 


z.lbo 


1 n 
l.U 


1 1 A O 

1.14Z 


Q 
O 


3 


o.Uoo 


7 


f.Uoo 


z.333 


o on a 


r 




A 07 


1.66 r 


1 con 


/I 

4 


r 




4.9b 


15 


14.9DO 


3.0 


o.Ulo 


1 n 
11) 


lU.Uyo 


z.O 


o no/? 
Z.UOD 


5 


o 

8 


o no r 

8.025 


on 

30 


on nor 

30.026 


3.75 


3.742 


on 

20 


1 n noo 

19.928 


2.5 


o ^ o o 

2.483 


« 
u 




1 9 QS^ 


uo 


^7 Q7Q 






oo 


OO.UOl 


9 Q9^ 


9 


7 


21 


21.01 


109 


109.015 


5.19 


5.189 


71 


70.964 


3.381 


3.378 


8 


34 


33.994 


201 


200.989 


5.912 


5.912 


130 


130.025 


3.824 


3.825 


9 


55 


55.004 


365 


365.008 


6.636 


6.636 


235 


234.983 


4.273 


4.272 


10 


89 


88.998 


655 


654.995 


7.36 


7.36 


420 


420.012 


4.719 


4.719 


11 


144 


144.001 


1164 


1164.004 


8.083 


8.083 


744 


743.992 


5.167 


5.167 


12 


233 


232.999 


2052 


2051.997 


8.807 


8.807 


1308 


1308.005 


5.614 


5.614 


13 


377 


377.001 


3593 


3593.002 


9.531 


9.53 


2285 


2284.997 


6.061 


6.061 


14 


610 


610.0 


6255 


6254.999 


10.254 


10.254 


3970 


3970.002 


6.508 


6.508 


15 


987 


987.0 


10835 


10835.001 


10.978 


10.978 


6865 


6864.999 


6.955 


6.955 



Table 3. Values of F n , T n , A n , C n>1 , and A n>1 for k = 3; <j> » 1.8392868 



n 


F n 


Appr. 


T n 


Appr. 


A n 


Appr. 


Cn,l 


Appr. 




Appr. 





1 


0.618 





0.132 


0.0 


0.213 





0.2 


0.0 


0.323 


1 


1 


1.137 


1 


0.946 


1.0 


0.832 


1 


0.75 


1.0 


0.66 


2 


2 


2.092 


3 


3.034 


1.5 


1.45 


2 


2.083 


1.0 


0.996 


3 


4 


3.848 


8 


7.96 


2.0 


2.069 


5 


5.126 


1.25 


1.332 


4 


7 


7.078 


19 


19.017 


2.714 


2.687 


12 


11.808 


1.714 


1.668 


5 


13 


13.018 


43 


43.028 


3.308 


3.305 


26 


26.095 


2.0 


2.005 


6 


24 


23.943 


94 


93.948 


3.917 


3.924 


56 


56.046 


2.333 


2.341 


7 


44 


44.038 


200 


200.032 


4.545 


4.542 


118 


117.891 


2.682 


2.677 


8 


81 


80.999 


418 


418.008 


5.16 


5.161 


244 


244.07 


3.012 


3.013 


9 


149 


148.98 


861 


860.968 


5.779 


5.779 


499 


499.007 


3.349 


3.349 


10 


274 


274.017 


1753 


1753.025 


6.398 


6.397 


1010 


1009.948 


3.686 


3.686 


11 


504 


503.996 


3536 


3535.998 


7.016 


7.016 


2027 


2027.043 


4.022 


4.022 


12 


927 


926.994 


7077 


7076.985 


7.634 


7.634 


4040 


4039.994 


4.358 


4.358 


13 


1705 


1705.007 


14071 


14071.015 


8.253 


8.253 


8004 


8003.979 


4.694 


4.694 


14 


3136 


3135.997 


27820 


27819.995 


8.871 


8.871 


15776 


15776.023 


5.031 


5.031 


15 


5768 


5767.998 


54736 


54735.994 


9.49 


9.49 


30956 


30955.993 


5.367 


5.367 



Since < l/4> < 1, this gives the asymptotic formula F n ~ <p n+1 /y/h — specifically, we have 



m+i 



where the error term 

E(n) 



1 /-1^ +1 
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Table 4. Values of F n , T n , A n , C n>1 , and A n>1 for k = 4; <f> « 1.92756198 



Tl 


F 

A n 




T 


A nnr 


A 


A nnr 


C 1 


A nnr 


A i 


A nnr 


o 




566 


n 


1 62 





287 


n 


0.212 





374 


i 




1 092 


x 


931 


1 

J. • W 


853 


1 


729 


1 


667 

\J . \)\J I 


2 


2 


2.104 


«j 


2 986 


1 5 


1.419 


2 


2.023 


1 


961 


O 


4 


A 056 


e 

O 


8 053 


9 


1 986 


<\ 

O 


5 091 


1 95 

.1. . _ "J 


1 955 


4 


8 

o 


7 81 8 

1 .010 


20 


1 9 951 


2 5 


2.552 


12 


12.11 


1 5 


1 549 


5 


15 


15.07 


47 


46.993 


3.133 


3.118 


28 


27.77 


1.867 


1.843 


6 


29 


29.049 


107 


107.033 


3.69 


3.685 


62 


62.063 


2.138 


2.136 


7 


56 


55.994 


238 


238.024 


4.25 


4.251 


136 


136.082 


2.429 


2.43 


8 


108 


107.931 


520 


519.932 


4.815 


4.817 


294 


294.017 


2.722 


2.724 


9 


208 


208.044 


1120 


1120.024 


5.385 


5.384 


628 


627.863 


3.019 


3.018 


10 


401 


401.017 


2386 


2386.03 


5.95 


5.95 


1328 


1328.068 


3.312 


3.312 


11 


773 


772.985 


5037 


5036.995 


6.516 


6.516 


2787 


2787.047 


3.605 


3.606 


12 


1490 


1489.977 


10553 


10552.957 


7.083 


7.083 


5810 


5809.98 


3.899 


3.899 


13 


2872 


2872.024 


21968 


21968.029 


7.649 


7.649 


12043 


12042.935 


4.193 


4.193 


14 


5536 


5536.004 


45480 


45480.014 


8.215 


8.215 


24840 


24840.053 


4.487 


4.487 


15 


10671 


10670.99 


93709 


93708.986 


8.782 


8.782 


51016 


51016.018 


4.781 


4.781 



approaches at an exponential rate. Indeed, this is the same as our asymptotic F n ~ 
4> n+1 /G'{a) for k = 2, because for k = 2 we have (p = (1 + V5)/2 and C7(x) = 1 + 2x, 
so G'(ct) = G'(l/0) = y/h. 

Finally, because all of our asymptotic formulas have absolute error terms that decay expo- 
nentially, in Tables [21 [3l and H] we obtain excellent approximations to our values F n , T n , A n , 
C n j, and A n j even for small values of n. 
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